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ABSTRACT 


Transmission  and  reception  cf  electromagnet ic  (Tv' 
energy  hy  communications,  weapons,  and  active/passive 
sensor  systems  is  known  to  he  strongly  influenced  "t y  an 
atmospheric  phenomena  known  as  ducting,  caused  ty 
refractive  layers  in  the  atmosphere  of  marine  environments. 
The  \'aval  Postgraduate  School  (NPS)  has  developed  a  Marin>- 
Atmospheric  Boundary  Layer  (MABL)  model  which  car.  he  used 
to  predict,  over  a  24  hour  period,  the  refractive  profile 
of  the  lower  atmosphere.  This  thesis  examines  the  co^el 
from  the  statistical/trend  analysis  approach  to  examine 
whether  the  model  can  he  used  as  a  valid  predictor  of 
refractive/duct  in g  conditions. 
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I . 

The  phenomena  mown  as  ducting  occur?  when  refractive 
layers  cause  electrorraenet ic  (FM)  energy  to  tend  toward  the 
earth  at  a  rate  greater  than  or  equal  to  the  curvature  o'* 
the  earth.  Ducts  occur  in  two  forms?  (1)  surf  ace-1  ased 
ducts  caused  ty  either  a  surf ace-ha sed  layer.  Figure  1,  or 
hv  an  elevated  layer,  Figure  2,  and  (?)  elevated  iu-ts 
caused  by  an  elevated  layer.  Figure  3.  As  shewn  in 
Figures  1  through  3,  duct  types  and  thicknesses  are  easily 
determined  ty  utilizing  the  profile  for  the  modified  index 
of  refraction  M  which  depends  on  temperature,  pressure,  and 
specific  humidity. 

The  Navy  employs  the  Integrated  Refractive  Fffects 
Prediction  System  (IRSPS)  to  identify  ducting  condition; 
and  to  assess  the  effects  on  various  fleet  emitter?  [Ref. 
l],  A  drawback  to  IF.EPS  usage,  from  a  tactical  standpoint, 
is  that  it  depicts  the  ducting  conditions  only  at  the  tire 
of  the  radiosonde  sounding. 

From  a  tacticians  viewpoint,  the  M-profiles  are  a 
valuable  output  of  the  Marine  Atmospheric  Boundary  Layer 
(MABL)  model  because  of  the  ease  with  which  duct  types  and 
thicknesses  are  obtained.  The  frequencies  which  are 
trapped  ty  ducts  ar®  a  function  of  duct  thickness,  as 
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igure  1 . 


Figure  2 
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Surface  Duct  Caused  by  Surface-Based  Layer 


Surface  Duct  Caused  by  Elevated  Layer 


A 


M  units 

Figure  3.  Elevated  Duct  Caused  by  Elevated  Layer 


stated  by  Eerr  (1951).  Thereby,  in-owing  the  duct  thickness 
and  type,  the  tactician  would  be  able  to  determine  the 
frequency  range  which  would  be  trapped.  In  order  for  the 
Ev  energy  to  be  trapped,  the  transmitting  antenna  must  be 
physically  located  within  the  duct.  Antenna  heights  are 
readily  obtainable  from  ship's  cha ra cte ri s t ic  cards,  hence 
the  transmitters  on  board  which  would  be  affected  by  the 
duct  could  be  quickly  identified.  The  ability  to 
accurately  predict  the  evolution  of  the  ducting  conditions 
for  a  24  hour  period  would  be  a  powerful  tool  in  the  hands 
of  today's  Carrier  Tattle  Croup  Commanders.  This 
information  could  be  used  to  establish  "mission  Control 
(IYCON)  policy  and  to  establish  effective,  effi-ient 
frequency  monitoring  plans.  The  M.A3L  Adel,  developed  by 
the  Environmental  Physics  Group  at  the  Naval  Postgraduate 
School  (NFS),  may  be  the  medium  through  which,  this 
information  can  be  made  available  to  the  fleet  tacticians. 

The  purpose  of  this  thesis  is  to  perform  a  detailed 
statistical/trend  analysis  of  the  MAPL  to  determine  its 
validity  in  predicting  ducting  conditions  for  a  24  hour 
period.  Cnee  this  analysis  has  teen  completed,  conclusions 
and  recommendations  will  be  made  in  order  to  provide  a 
foundation  for  further  refining  of  the  VA3L  model,  if 
necessary,  by  the  Environmental  Physics  Group.  The 


overall  objective  is  to  prcvile  a  reliable  predictive 
model ,  available  for  dissemination  and  incorporation  in  * ve 
fleet  . 

The  approach  taken  in  this  study  differs  froT  pest 
MABL  rodel  evaluations  at  NFS,  in  that  this  approach  Is  t c 
perform  an  objective  examination  of  the  Todel .  The  goal  is 
to  extract  information  about  the  model  output  from  a  purely 
statist  ical/trend  analysis  approach,  which  has  net  re-'’ 
done  previously. 
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II.  *OD£L  iISCRI£2IO\ 


Ths  ''l  A  PI  rr  od  e  1  ic  a  ^“rc— order,  two  1  a  y*  p  r »  L  z  *  ~  A  *  3  c1 
mixed  layer  model.  The  atmosphere  is  assumed  to  corsist  of 
two  layers:  a  well-mixed,  turbulent  boundary  layer,  and  the 
relatively  non-turbulent  free  atmosphere  above.  The 
inversion  ( cr  "transition  zone”)  separates  one  layer  frcT 
the  other.  In  a  zero- order  model,  the  inversion  thickness 
is  assumed  to  be  zero,  hence,  a  discontinuity  or  “juri" 
occurs  in  the  profiles  of  the  variables  at  the  inversion 
[ref.  ?]. 

The  NFS  model  of  the  rrixed-layer  is  implemented  on  a 
Fevlett  Packard  9845  desktop  computer.  The  model  uses  a 
32-ninute  tire  step  to  predict  mixed  layer  temperature, 
specific  humidity,  the  jump  of  these  values  at  the 
inversion,  and  cloud  or  fog  formation. 

The  M  profile  is  calculated  from  predicted  r  and  : 
values  using  Equation  (1). 

5  2 

M  =  7?.S(P/T)  +  6  x  10  *  (q  *  P/T  )  +  2.157  *  ?.  fi' 

where  ?  is  pressure  in  millibars,  T  is  temperature  in 
degrees  Kelvin,  q  is  specific  humidity,  and  Z  is  height  in 
meters.  The  top  of  the  duct  corresponds  to  the  height 
above  the  surface  wh^re  the  y  value  is  minimum.  The  base 


it-  first  intersects  a  point  of  equal  v-units  or  the 
surface,  whichever  occurs  first.  ducting  ~ err  only  =  r ;  = ' ; 
where  there  is  a  ter pe mature  inversion  which  acrs  a;  a 
trapping  layer,  refracting  ?M.  energy  toward  the  earth. 
Inversions  typically  exist  in  rarine  surface  high  pressure 
regions  where  well-mixed,  stable  layers  are  found  between, 
warm  dry  air  above  and  cooler  mist  air  below.  'a"?»  •-**£ 
of  low  level  stratus  clouds  often  indicate  areas  of  \  i;  :•  t 
occurence.  Ducting  is  expected  to  be  minimal  n-ar  f-ents 
and  areas  of  convective  cloudiness.  -rents,  with  their 
associated  upward  motion,  often  dissipate  the  inversions  as 
the  whole  air  column  becomes  mixed.  4reau  of  convert iv- 
activity,  discernible  by  the  presence  of  cumulus  clcds, 
are  also  normally  inversion  free  [P.e f.  1]  . 


III.  C^IJEF_IA  ANE  ASfUMfT  ICN£.  F££  MABL  ”£EEL  V EDIFICATION 

A  major  assumption  in  the  verification  of  this  ~ooei 
is  that  the  radiosonde  data  is  correct.  This  point  will  he 
elaborated  upon  later  in  this  paper. 

A.  CRITERIA  FOR  CHCOSINC  EXAMINATION  PIRIOES 

Forty  cases  (ex ami  nation  periods'  were  chosen  to 
evaluate  the  V,AEL  model.  The  radiosonde  soundings  use:  in 
all  4?  cases  were  tak’en  off  the  coast  of  Foutherr. 
California  in  end  around  fan  Clemente  Island.  These 
examination  periods,  which  are  2 4  hours  in  length,  were 
cho^p  based  on  the  following  criteria: 

1.  There  had  to  be  at  least  one  radiosonde  sounding 
ir.  the  25  hour  period  preceding  t^e  examination  period  with 
surface  layer  data,  from  the  NFS  research  vessel,  Alania, 
within  ?.  hours  of  that  sounding.  This  requirement  is 
needed  to  calculate  Vs  (subsidence)  with  Lenschcw's  (137o) 
method . 

2.  There  had  to  be  surface  layer  data  within  ?  hours 
of  the  radiosonde  sounding  at  the  beginning  of  the 


examination  period,  also  required  in  calculating  *s . 


i  r.  a  d  d  i  t,  i  o 


3.  There  had  to  he  at  least  z  s ound i ngs , 
to  the  initial  sounding,  during  the  examination  period  to 
he  used  as  verifications  of  the  model's  predictions. 

4.  Wind  speeds  had  to  he  available  within  3  hours  of 
the  initial  sounding.  If  wind  speeds  are  available  within 
3  hours  of  the  12  hour  point  or  within  3  hours  of  the  24 
hour  point  they  will  also  he  used.  If  these  2  values  *  re 
not  available,  the  initial  value  will  he  used  and  held 
constant  during  the  exam  period. 

5.  INITIALIZATION  RETIREMENTS  FOR  yOT  EL  P  -  ELI  CT  10 \S 

The  model  requires  as  input  data  the  following? 

1.  Radiosonde  Sfiuiii Lg  Ldii 

Paw  data  from  the  radiosonde  launches  must  be 
digitized  for  use  in  the  VAEL  model. 

2.  Suh s iden eg  {.W«l  Prgdigtion 

The  subsidence  value  used  in  the  model  prediction 
for  all  4 2  cases  was  obtained  from  integration  of  tv> 
Moisture  Budget  Equation  given  in  Equation  (2).  This 
method  was  developed  and  performed  by  Lenschow  (197?).  it 
is  based  on  the  assumption  of  well-nixed  specific  humidity 
(q)  in  the  boundary  layer,  and  that  changes  occur  due  to 
fluxes  at  the  sea  surface  and  inversion  only. 

**  *  ala* i)  »h/<n)  -  hUic*n/ni  + 


rA\W*l 


liquid  water  contents,  w'q'  is  the  moisture  flux,  /3  is  t hc 
vert  ice  i  gradient  of  (~-»l;  in  the  mixed  ’a  ;‘-r,  non  A-'c-i) 


is  the  difference  between  total  moisture  in  the  mixed  layer 
and  total  moisture  immediately  above  the  ,mixci  lay^v.  In 
this  evaluation,  the  linear  trend  in  subsidence  during  the 
24  hour  period  prior  to  an  examination  period  was 
calculated,  and  the  regression  line  applied  at  the  model 
start  tire  to  find  the  suhsii°nce  value  used  th roughen*  the 
prediction  period. 

If  the  subsidence  value  was  positive.  then  a  d0 fault 
value  of  zero  was  used  to  initialize  the  model.  A  positive 
value  of  subsidence  causes  the  model  to  "hlow-u^".  The 
reason  for  this  has  not  been  determined  at  the  tire  of  this 
writing. 

7.  Tea  Supfcge  Te2£erat.ure  LZ5.21 

This  value  was  obtained  from  Acania  data.  The 
SST  at  the  beginning  of  the  period  was  held  constant 
throughout  the  examination  period. 

4 .  Wind  Speed 

a.  If  wind  speed  was  not  available  within  7  hours 
of  the  12  hour  point  but  is  available  at  24  hour  point, 
then  a  linear  regressi  on  was  performed  using  the  i n  i 1. 1  a  1 


and  24  hour  values 


17.  STATIST  I CAl/lRENT  AN ALIBIS  07  THJ  ^ A7L  *2111 

A.s  previously  stated,  the  main  objective  of  this 
thesis  is  to  determine  the  valid! ty/accuracy  of  the  VA?L 
model's  ability  tc  predict  ducting  conditions  ever  a  ?A 
hour  period  . 

Of  the  42  cases  used  in  the  verification,  five  ''ases 
were  omitted  because  during  the  examination  periedc  the 
atmosphere  was  stable,  therefore  the  model  physics  dc  "ot 
apply.  Seven  cases  were  omitted  because  t^e  inversion 
height,  Zi,  was  less  than  222  meters,  for  which  ^ase,  tie 
model  physics  do  not  apply.  This  left  28  cases  to  be  used 
for  the  model  verification .  These  28  cases  contain  a  lota] 
of  89  radiosonde  soundings  used  for  verification  purposes. 
Although  the  model  predicts  several  parameters 
(T,q,Zlcl,Zi  )  during  the  24  hour  prediction  period,  the 
parameter  under  evaluation  in  this  thesis  is  th«- 
temperature  inversion  height,  Zi ,  because  this  value  is  the 
height  of  the  top  of  the  atmospheric  duct.  If  this  value 
cannot  be  predicted  accurately,  the  model  is  of  little  use 
to  the  fleet  tactician. 

The  current  method  employed  by  C-eouhysics  rffictr? 
afloat  to  predict  T  (temperature),  q  (specific  humidity), 
Zlcl  (lifting  condensation  level),  and  Zi ,  is 


"Persistence”.  Persistence  means  that  nothing  changes. 
Per  example,  if  a  radiosonde  launch  at  <?S.d?  indicates  that 
the  inversion  height  is  meters,  then  a  persistence 
prediction  would  state  that  the  inversion  height  would 
remain  unchanged,  until  the  next  launch  at  which  tire  it 
would  he  updated  to  the  current  value. 


4.  STATIST I  CAL /TP EM  ANALYSIS 

TATA  AND 

-PS’JLTS 

The  procedures  and 

methods 

use"  for 

*  hr 

s ta t i s t ica 1 /trend  analysis 

of  the 

YABL  red  e  1 

s  n ' 

persistence  will  now  he  presented..  All  of  tve  results 
presented  in  this  chapter  will  he  discussed  in  the  chapter 
of  conclusions. 


1  •  rue t  ll£ nd  Pred ic t ip:} 

Prior  to  examining  the  inversion  height 
prediction,  it  is  important  to  determine  how  accurately  the 
'lABL  model  and  persistence  predicts  the  trend  of  existing 
ducts  during  the  24  hour  examination  period.  Predicting 
the  trend  hy  the  model  means,  for  example,  if  an  elevated 
duct  occurs  at  the  initial  sounding  a  no  tlev  tsoTcs  a 
surface  based  duct,  does  the  model  predict  this  tr^udL 
Predicting  the  trend  hy  persistence  involves  looking  at  thc 
radiosonde  launch  immediately  prior  to  the  initial  sounding 
and  observing  the  trend.  For  example,  if  the  inversion 
height  is  decreasing  during  this  period  then  it  is  assured, 


IS 


continue 


decrees® 


4 


exa~i  net  ior.  period.  This  segment  of  the  evaluation 


not  consider  the  accuracy,  in  meters,  of  the  predict 
only  the  trend.  This  means  the  r  ®  could  be  a  si.;r.ifi 
differ®nce  at  the  verification  time  between  Zi  of 
radiosonde  a"d  Zi  of  the  model  or  persistence,  yet 
trend  could  hay®  teen  correctly  predicted. 

Considering  all  28  cases  un^er  examination, 
"•cdel  predicted  the  trend  correctly  "Z.^X  of  the  t 
whereas,  persistence  predicted  the  trend  correctly  in 
of  the  <*a 5 es  . 


Th®  25  cases  were  further  subdivided  i^to  c 
or  cloudy  sky  conditions,  existing  at  the  onset  o* 
examination  period,  to  if  the  initial 

conditions  had  any  affect  on  the  trend  predict!  or.s.  T 
were  IT  cloudy  and  1?  clear  sky  cases.  Both  the  mcdel 
persistence  predicted  the  duct  trend  correctly  in  7ri  ( 
1?)  of  the  clear  sky  cases.  Persistence  was  correct  in 
(?  of  lc)  of  the  cloudy  sky  cases,  however,  the  nod®! 
correct  in  55. (12  0'  IS)  of  th®  cloudy  s.-cy  ces®s. 


data  is  shown  in  Table  I 


I 


TAr'II 

Duct  Trend  Prediction 


•r  Vie  p 
•-J  -t  —  i- 

hk* 

^2l  fl 

TCorrec t 

Persist  ^nc6- 

*  2  j  r 

25 

A  1 1 

17  correct 

50.7 

If  correct 

it  <~l 
_  r 

1  c 

^  ^  1 1  A  \ r 

*  -  W  •  •*  C 

12  correct 

56.7 

?  correct 

e 

1 2 

Clear 

7  correct 

70.0 

7  correct 

70 

2.  Inversion  Height  Lli.1  ?L££i  CilSIl 

Tatle  III  is  a  complete  listing  of  the  raw  date 
sed  in  the  analysis.  This  taMe  gives  the  Tate,  C^s  = 
urber,  ntirrter  of  hours  since  the  initial  radiosonde 
cunding,  the  radisonde  inversion  height,  VA?L  rrodfl 
redictlon  of  Zi,  and  persistence  prediction  of  Zi. 
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Conplete 

Data  Listing 

for  Z i 

C  5.5  <?£ 

“aurs. 

Eiii os  end  e 

V  p  1 

.  'j  •*.  r* 

Persist0 

Initial 

25  Sep. 

2 

12 

1445 

1 559 

1 455 

1975 

25 

135* 

1497 

1446 

25  Sep. 

'Z 

i,. 

11 

12°4 

1723 

1265 

1975 

24 

114 

2036 

1  294 

03  Oct  . 

7 

3 

327 

224 

47^ 

197« 

21 

154 

54 

307 

24 

150 

65 

154 

10  Cct. 

10 

11.5 

193 

192 

250 

1975 

24 

336 

142 

193 

11  Oct. 

11 

11.5 

237 

555 

'Z  7<? 

1975 

22.5 

455 

770 

237 

20  July 

12 

15.5 

ise 

747 

510 

1977 

24 

305 

671 

195 

21  July 

13 

l*. 5 

190 

4  ?  P 

•*  ?  ' 

TABLE  III  (c oat.) 


Tate 

Cas  e* 

Eiii2£0Ri_e 

:L2£el 

Pc  r> c  i  c 

ifci.  S»£.  ** 

Initial 

22  1^3  v 

24 

4 

922 

732 

5tc 

1975 

12.5 

733 

=  25 

372 

23.5 

1515 

955 

7?  3 

71  Aug. 

25 

14.5 

352 

16  = 

475 

1975 

2%  c 

3  =  4 

112 

729 

w  *— 

22  Km. 

2  £ 

12 

429 

424 

197  = 

24.5 

352 

461 

479 

?3  Aug . 

27 

12 

31  = 

337 

^  .  «. 

1072 

23 . 5 

352 

345 

31  = 

24  Aug  . 

25 

12.5 

2  99 

52  7 

35? 

1973 

24 

241 

F  1  2 

35  Au?. 

29 

12 

197 

235 

Pd  1 

1979 

24.5 

165 

232 

197 

17  9  Aug. 

32 

12.5 

319 

314 

23? 

197£ 

19 

354 

335 

315 

24.5 

334 

357 

364 

23  Aug. 

-  <»- 

C 

TTC 

-  w 

255 

734 

1975 

13 

7/7 
c  -c 

254 

77“ 

w  >». 

15.5 

445 

234 

343 

24 

454 

222 

44=’ 

12  Aug. 

34 

F 

326 

242 

314 

1379 

14 

233 

231 

24 

322 

225 

233 

11  Aug. 

35 

12 

555 

3S5 

322 

1973 

24 

S53 

435 

556 

12  Aug. 

*te 

\j  - 

12 

2  =  2 

1297 

=45 

19'7~ 

24 

=  57 

1435 

2  =  7 

13  Aug. 

37 

12 

31.3 

365 

cc.7 

197  = 

15 

215 

1237 

317 

24.5 

135 

1233 

216 

/*  .  '  '  m  y  V-V-  '-V-  V*  .  •  ‘  *  •-  /.*.  V  ”  v\  ‘\*  V%‘  ,V"  -V  ’; 


Hour«  Since  ^odel  ?ers: 

Initial 


Aug. 

1  c 

F 

792 

254 

09  z 

?c 

12 

512 

25  2 

A 

75? 

345 

F 

^  X 

May 

i2 

7 

45  2 

532 

54 

se 

15 

252 

73  = 

/  c 

r>T  c= 

525 

£  Fn 

pc 

—y 

43 

%  5 

352 

255 

0  0 

A*-*  X. 

52 

1? 

CR" 

7  ; 

2: 

e 

21 9 

Mav 

44 

c  e 

•  ■- 

954 

c.C'? 

c7‘ 

32* 

24.5 

254 

527 

9r 

In  order  to  perform  statistical  analyse 
the  predictive  accuracy  of  the  VA3L  model  and  persis 
the  tTewlett-Pacicard  HP-S5  was  used  tc  obtain  nec 
statistical  values.  The  prortre-’-  "PIST?"  was  vs 
assemble  inpu  +  data  sets  in  ascending  crde^  a' 
calculate  the  mean,  standard  deviation,  and  the  ra 
the  data  set.  The  program  PIN  was  used  to  divide 
data  set  into  bios  of  equal  width  to  form  a  histowr 
the  data  set.  "T°LCT"  was  then  used  r0  plot  each  dat 
The  program  PLOTP”  was  used  to  enter  data  pairs, 
"TPLCT"  was  used  to  make  a  scatter  plot  of  the  = 
points.  These  programs  are  listed  in  Appendix  a. 

In  an  effort  to  obtain  more  useful  inform 

?b.p 


five  particular  situations  were  evaluated 


~stercries  chosen  «er e:  l)  all  ?3  date  points,  2'  o’7  lata 
points  with  the  h.i^h  and  low  values  frr  category  1 
ercluded,  ’ )  the  ti^e  when  the  verification  was  taben 
during  the  2-  hour  examination  period,  V:  the  season  of  th ‘ 
year  when  the  verificatcn  was  taVoen,  and  c)  the  initial  shy 
conditions.  Within  categories  ?,  4,  and  f  there  were 

?.  ,  and  2  classes,  respectively.  Classes  within  os  tc..;  ori  r  s 
are  rutually  “^elusive.  Categories  1  and  7  had  only  one 
class  each.  In  all,  there  v e  r  e  clevc,>.  classes  e  v  a  1 1  ■  a  *  \  ~  , 
each  class  c on t  a i n inf  two  data  sets,  one  for  rods',  data  and 
one  for  persistence  data.  The  c  categories  and 

associated  classes  are  described  in  Table  IV. 

The  Tel t a  M  and  Delta  P  values  for  each  data  set 
are  listed  in  Appendix  2,  in  the  output  ferret  cf  the  -’?-ct 
"DISTF"  program.  Delta  Y  and  Delta  ?  were  calculate'5  vy 
subtracting  the  roiel/pers  ist  e.nee  prediction  frc'r  the 
radiosonde  value  at  that  tire.  The  histograms  cf  the 
prediction  errors  are  in  Appendix  C. 

Table  7  is  a  corpilation  of  the  data  set 
filenames,  the  r.unter  of  data  points,  N,  in  the  set,  an"1 
the  mean  and  standard  deviation  of  the  set  in  peters. 

Ideally,  the  desired  error  between  the 
radiosonde  sounding  and  the  prediction  value  (redei  cr 
persistence)  is  zero.  So,  the  desired  mean  o'  each  data 
set  is  zero.  The  desired  standard  deviation,  of  the 
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prediction  error,  is,  ideally,  plus  or  minus  5£  reters 
(plus  or  minus  100  meters  will  also  be  examinee' .  Any 
inversion  heights  which  cannot  be  predicted  within  plus  or 
minus  50  meters  is  of  little  or  no  value  to  c  tactician. 
Thi's  predicted  value  is  critical  because  frequencies  which 
are  trapped  within  a  duct  are  a  function  of  duct  thickness 
(Kerr  1951)  and  the  trao srrit t ing/r ecei ving  antenna  must  be 
located  within  the  duct  for  maximum  exploitation. 
Therefore,  the  standard  deviation  is  an  important  statist!-" 
in  that  it  is  a  measure  of  the  precision  of  the  prediction. 
In  other  words,  it  reveals  hew  much  fluctuation,  about  the 
mean,  occurs  for  a  given  data  set.  The  larger  the  standard 
deviation,  the  less  precise  the  prediction. 

The  given  mean  and  standard  deviations  do  not 
clearly  reveal  the  predictive  worth  of  either  the  mod*!  or 
persistence.  In  order  to  make  more  concise  statements 
about  the  two  predictive  methods,  it  will  be  useful  tc 
perform  hypotheses  tests  utilizing  this  data  and  a  level  of 
s igr.i f i cance ,  a  ,  of  0.01. 

a.  Hypothesis  Concerning  Cne  Pean 
Hypothesis  testing  is  done  to  determine  whether  or 
not  a  given  statement  concerning  a  data  set  can  be  accepter, 
with  a  probability  a  of  risking  a  Type  I  error  [Ref.  3:  p. 
195].  This  is  referred  to  as  the  null  hypothesis,  Ho. 


TABLE 

IV 

Category/C  lass 

Description 

Cs.ieg.crv  busier 

Class  Nurrter  Description 

1 

1 

All  53  data  points  use 

2 

2 

57  data  points  uscd. 
The  high  and  lew 
values  frorr  data  set  1 
were  excluded. 

T 

V. 

Badiosonde  soundings 
which  occurred  4  to 
12.33  hours  after 
initial  launch 

3 

4 

Badiosonde  soundings 
which  occurred  11  tc 

13  hours  after 
initial  launch 

3 

5 

Badiosonde  soundi^**? 
which  occurred 
to  22.39  hours  after 
initial  launch 

3 

Badiosonde  soundings 
which  occurred  23  to  2: 
hours  after  initial 
launch 

4 

7 

Radiosonde  soundings 
taken  during  ^arch, 
April,  and  May 

4 

Q 

Radiosonde  soundings 
taken  during  June, 

July,  anc  August 

4 

O 

Radiosonde  soundings 
taken  during  Septercer, 
October,  and  hcverter 

5 

13 

Sky  conditions  clear 

e, 

11 

Sky  conditions  cloudy 

and  Standard  re viatica 


veaa 


ate  Set 

N 

Mean 

Standard  Pe 

ilenare 

vodc 1  i 

53 

82. £12 

4 12 . 226 

Persistl 

£*  -V 

C  3 

~4 .565 

r-  O  /■“'  rj  C  ^ 

c  .  7Cc 

^od  e 12 

6? 

62.826 

335 .425 

Persist2 

67 

-11 .881 

217 .242 

vodel2 

12 

-121.122 

138.291 

So rc i c  t? 

12 

-57.622 

183.94? 

Model4 

15 

149 . 256 

.  251.582 

Pars  i  s  t4  ' 

15 

93.275 

•>  ?  z  -  " 

J 

y  0  C  9  1  £ 

12 

158.917 

*7-1-1  *7*777 

ecu#  f  V. 

Persists 

12 

42.917 

145  .  =  ?£ 

YodelS 

2? 

32.724 

363.166 

Persists 

27 

-88.519 

223  .396 

vode 17 

22 

55 . 554 

23? . 2®2 

Persist? 

22 

-48.591 

223  .03? 

Moc!elS 

35 

51.285 

392 .758 

Persist® 

35 

0.020 

241 .952 

Model? 

12 

81.922 

193 .225 

Persist? 

12 

27.320 

124.96? 

vode 110 

25 

22.922 

457.560 

Persist  1? 

25 

-36.720 

335.522 

'lodelll 

44 

142.223 

372.330 

?e  r  s  i  s  1 1 1 

44 

7 .870 

246.600 

REPRODUCED  AT  DO  V  E  R  N  M  E.N 


Figure  4.  Type  I  Frrcr 

Given  the  distribution  ir.  Figure  4,  a  Type  I 
error  would  be  the  probability  a  of  rejecting  the 
hypothesis  when  it  is  true. 

The  null  hypothesis  chosen  was.  Ho:  n  -  mc, 
which  in  words  states,  is  the  mean  of  Delta  '■'/Delta  ?  equal 
to  mo,  a  specified  value  for  the  desired  mean.  Since 
D-plta  m  and  Delta  P  are  the  errors  associated  with  the 
difference  between  the  model  or  persistence  prediction  and 
the  radiosonde  value  and  the  desired  error  to  be  zero,  then 
Mo  was  chosen  to  equal  zero.  Another  way  of  stating 
this  would  be:  Can  it  be  said  with  reasonable  confidence 
(99?)  that  the  mean  of  the  error  for  each  predictive  method 
be  assumed  to  be  zero? 

In  order  to  carry  out  this  test,  the  Student's 
t-test  was  used.  The  assumptions  required  for  this  test 
are  that  the  actual  standard  deviation,  o  ,  be  unknown  and 
that  the  sample  comes  from  a  normal  population.  The  t 


statistic  is 


REPRODUCED  AT  GOVERNMENT  EXPENSE 


X  -  4  0 

*^T 

where  x  equals  the  mean  of  the  sanple,  (jo=  £.0,  S  is  the 
sample  standard  deviation,  and  n  is  the  number  of  data 
points  in  the  set.  This  calculated  value  of  t  was  compared 
to  a  value  obtained  from  a  t-statistic  table.  This  value 
was  obtained  by  entering  the  table  with  the  number  of 
degrees  of  freedom,  in  this  case  n-1,  and  a/ 2  -  0.2Z5 
because  this  is  a  two-tail  test.  If  the  absolute  value  of 
t-calculated  is  greater  than  t-table  than  the  hypothesis  is 
accepted.  Figure  5  shows  the  acceptance/rejection  regions 
for  a  two-tail  test  [Ref.  3:  pp.  211-214]. 

The  results  of  the  hypothesis  test,  Ho:  p  =  no  for 
each  data  set  are  shown  in  Table  VI. 


Figure  5.  Acceptance/Rejection  Regions 


ABLE  VI 


F  o :  M  =  Mo  =  2 .2 


Eilenare 

n 

izsalc 

iriftilfi 

£-£S£it  ZEf  Jgoi 

M odell 

59 

1 .640 

2.5^6 

A  cceot 

Persistl 

69 

-2.135 

-2.5*5 

Accept 

wocel? 

57 

1  .533 

2.576 

Accept 

Pers is  t2 

57 

-0.44S 

-2.576 

Accept 

*  od  e  1? 

1? 

-1  .930 

-^.250 

A  C  n  6s  J' t 

Persist? 

1? 

-0.990 

-3.250 

Accept 

vodel* 

18 

2.415 

2.698 

Accent 

Per  s i s  t  4 

1 8 

1 .688 

2.898 

Accept 

vodel5 

12 

1  .650 

3.105 

Accept 

Persists 

12 

0.952 

3.106 

Accept 

^odelS 

27 

0.416 

2.779 

Accept 

Dersist5 

27 

-2.05? 

-2.779 

Accept 

yodel7 

22 

0  .926 

2.611 

Accept 

Persist? 

22 

-1.123 

-2.531 

Accept 

vodelS 

35 

0  .927 

2.576 

Accept 

Persists 

35 

2.000 

2.576 

Accept 

Model9 

10 

1  .307 

3.250 

Accept 

Persist? 

10 

0.691 

3 . 2  50 

Accept 

Mode  115? 

25 

0  .254 

2.797 

Accept 

Persistl? 

25 

-2  .577 

-2. 797 

Accept 

Modelll 

44 

2.490 

2.576 

Accept 

Pers istll 

44 

0.210 

2 . 575 

i  mpTit 

t .  Hypothesis  Concerning  One  Variance 
The  null  hypothesis  considered  in  this  section 
will  he  rc  determine  if  the  population  standard  deviation 
is  equal  tc  a  specified  constant  (Ho:  o-  a  o).  This  test 
is  conduct  ai  ty  cal  evict  in.?  the  Chi-Square  statistic, 


u-t>s2 

o\ 


where  n  is  the  number  of  data  roints  in  the  sarnnle 
2  ‘  2 
population,  ?  is  the  variance  of  the  sarple,  and  err  is 

the  specified  (desired)  variance  (variance  is  the  square  of 

the  standard  deviation).  The  Chi-Square  calculated  value 


is  then  corps  red  to  a  Chi-Square  value  obtained  from  a 
table  in  Ref.  3,  using  n-1  degrees  of  freedom  and  a~  ?.P1. 
If  the  calculated  value  is  less  than  Va/^  or  greater  than 
null  hypothesis  is  accepted.  The  desired 
standard  deviation  was  chosen  to  be  °  o  =  ~,Z  meters, 

however,  °o  =  12(?  neters  was  also  tested.  The  results  of 


these  tests  are  given  in  Table  V'll  [Pef.  3:  op.  233-235]. 

c.  Hypothesis  Concerning  Two  Means 

In  order  to  determine  if  there  is  a  significant 

difference  between  the  two  means,  the  model  and 

persistence,  the  null  hypothesis,  Ho:  pm  =  pp  must  be 

examined.  Along  with  the  assumption  of  normality,  a 

statement  must  be  made  about  the  variances  of  the  two  data 

2 

sets.  Namely,  can  it  be  assumed  that  the  variances  am 
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[IV  -V  ■  *V  ’T,'  "■  ■n.  ■v 


”WV  'V^T' 


TAEL2  VII 


Ho 

:  o  = 

o  o 

ao  =  5? 

and  100 

meters 

2 

2 

X 

?i  1 e  o  a  me 

S 

X’ziiiS 

ulaiai 

X=i  iii£ 

hZ2&ii 

/o  £  i  O  f  t 

K 

122 

50 

100 

Moiell 

1.7X1?“ 

4. 

4570 

1140 

102.3 

P  P  iPCt 

Re  'ect 

Persistl 

7.9x10 

c 

213  0 

532 

102.3 

Reject 

~  n  i  £>  r*  t 

V-  ^  W  ^ 

v  od  9  1 2 

1.1x10 

4 

2980 

746 

99. 5 

Reject 

^  r  t 

Persists 

rr 

4.7x10 

£ 

1250 

312 

cc  a 

•  /  w  •  V 

E  e  .5  e  c  t 

?  e  j  e  c  t 

M  cde 13 

3.9x10* 

141 

T  a 
-• 

23.6 

Reject 

pL0  i  £  ~  t. 

Persi st3 

3.4x10 

122 

30 

23.5 

F  e . j  e  c  t 

"e  if"  ; 

Model4 

9.9x10 

466 

117 

35.7 

Reject 

Fej  ec  t 

4 

Persi s  t4 

4.4x10 

e. 

299 

?4 

35.7 

Reject 

Reject 

M  od  e  1 5 

1.1x10“ 

488 

122 

26.8 

Reject 

?.  e  '  £  ■>  t 

Persists 

it 

2.2x10 

e 

97 

24 

25 .  S 

Fej  ect 

Ac  cep  t 

vodelS 

1.5x10 

i 

1530 

3S2 

48.3 

Reject 

Fe jec  t 

Persist? 

5.0X10* 

4 

522 

131 

46.3 

Reject 

F.e  jec  t 

vodel7 

8.3x10 

4 

693 

1?3 

41.4 

F.e  j  ec  t 

Reject 

Persist? 

4.1x10 

By 

345 

96 

41.4 

Re jec  t 

Re jec  t 

Models 

1.5x10 

4 

2070 

517 

56.1 

Reject 

Re jec  t 

Persists 

5.9x10 

4 

796 

199 

56.  1 

Fej  ec  t 

Fe jec  t 

Model9 

3.9x10 

141 

35 

23.6 

Reject 

Reject 

4 

Persi st9 

1 .9x10 

5 

55 

14 

23.6 

Reject 

Accept 

Mode  110 

2.0x10 

1950 

487 

45.6 

Re  jec  t 

Re jec  t 
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TABLE  VII  (cont .  ) 


Filename 

a 

2 

X-oalcula£ed 

~W~  ill 

2 

Xriabls 

Accent /Reject 
~-50~"l30 

Persistl0 

- / 

1.1x10 

10S3  271 

45.5 

F.e  jec  t 

F.e  jec  t 

■'•lode  111 

1.4x10 

2?9?  559 

71.0 

Re jec  t 

P.e  jec  t 

Persistll 

6.2x13 

1060  265 

71.0 

Re.iec  t 

Fe jec  t 

2  2  £ 
and  <*z  are  equal?  The  null  hypothesis  Ho:  =  ap 

must  he  completed  to  answer  this.  If  it  can  be  said  that 

they  are  equal,  (i.e.  accept  the  hypothesis),  then  a  two- 

sample  t  test  can  he  used  to  test  the  means.  If  the 

hypothesis  that  the  variances  are  equal  is  rejected,  then  a 

paired-sample  t  test  must  he  used  to  test  the  means. 

In  order  to  perform  the  hypothesis  test,  an  F- 
statistlc  was  calculated  hy. 


2.  a 

where  Sm  is  the  variance  of  the  model  data  and  Sp  is  the 

variance  of  the  persistence  data.  Using  a-  0.01,  if  F- 

A 

calc  was  less  than  J-table  then  the  hypothesis.  Ho:  or.  - 
ap  must  be  accepted.  The  results  are  shown  in  Table 


-yji  n-1 )  [Sn  -- 

a  =  ? .  ?i  and  the  decrees  c 
esulis  are  giver,  in  Tetle  I) 


T  r r  L I  I> 


1.2x1? 


5. 5x 1? 


3.2x10 


1.37 

0.70 

0.22 


In  order  to  conduct  this  test 
values.  Delta  M  -  Delta  P. 


2.5? 

Accept 

2.59 

Ur-'O'  t 
*  -  -  -  ±*  v 

2.32 

Accept 

2.59 

Accept 

for  not 

being  ahl 

~es  are 

equal  had 

e  paired 

-sample  t 

data  set 

contain i n  g 

compiled 

and  a  mean 

sets 


as  in  C°cticn  A 


The  six  paired 


sets 


a  r  e 


listed 

in 

Appendix  P  ard 

the  results  of  the 

hypotnesi s . 

Ho:  fx 

= 

H o  =  0.0,  are 

listed  in 

Table  X. 

taplf  x 

Fc: 

n  =  no  = 

0.0 

Pa  t  a 

^  o  t 

£  Izc^ipula 

ted  iz 

table 

Acre  i;  t  1  ?.  e  4 

1 

2.31 

• 

2  .££ 

Reject 

2 

3.01 

2.5  = 

Reject 

1  .£? 

3.06 

Ac  cep t 

*0 

2.53 

2.75 

Accept 

- 

1 .  £4 

Accept 

ii 

3.73 

p  xp 

Reject 

d.  Prediction  Error  vs.  Verification  Time 
Figures  £  and  7  are  plots  of  the  Prediction  Frror 
vs.  Verification  Time.  :*'ith  zero  being  the  desired  error, 
it  is  clear  that  as  the  model  progresses  further  into  the 
?4  hour  prediction  period,  that  the  prediction  error 
increases  significantly.  The  standard  deviation  varies 
from  13?. 031  meters  to  3-7 .155  meters  (from  Table  V). 


VERIFICATION  TIME  (HOURS) 


PERSISTENCE  PREDICTION  ERROR  VS.  VERIFICATION  TIME 


VERIFICATION  TIME  (HOURS) 

Figure  7.  Persistence  Prediction  Error  vs.  Verification  Time 


The  Persistence  Prediction  Freer  plot  reveal 
much  less  scattering  of  the  data  points,  while  trer?  is  n 
obvious  trend  to  the  prediction  error. 

e .  Pr ed  i c  ti  on  Fr r  or  vs  .  - c n t h  of  Ob  se  rve  t  i  or 
Figures  8  and  5  are  plots  of  the  ??  data  point 
relative  to  the  month  in  which  they  ere  taken.  The  y-axi 
of  these  two  firures  are  labeled  by  numbers,  where 
corresponds  to  January  through  12  which  is  December.  It  i 
obvious  that  there  is  much  more  scatter  associated  with  tv 
model's  prediction  error  during  the  summer  months  then  = 


other  times. 


MODEL  PREDICTION  ERROR  VS.  MONTH 


Prediction  Error  vs.  Month  of  Observat 


stence  Prediction  Error  vs.  Month  of  Observat 


V.  CONCLUSIONS 


A.  DUCT  TREND  PREDICTION 

Erorn  a  statistical  viewpoint  there  is  no  significant 
difference  between  the  MABL  model  and  persistence  in 
predicting  the  duct  trend  or  sky  (clouds/no  clours) 
conditions.  There  are  seme  differences  in  their  correct 
prediction  percentages,  however,  net  enough  to  cheese  one 
method  ever  the  other.  The  results  clearly  indicate  that, 
given  sky  conditions  are  clear  at  the  initial  start  time, 
the  duct  trend  car.  be  accurately  predicted  by  either 
method.  Cloud  formation  appears  to  c^^^pr cblems  in  the 
model's  predictive  ability  and  should  be  further 
i nvest igat ed . 


E .  HYPOTHESES  TESTS 

Although  the  hypothesis  test  of  Section  Ilia,  showed 
that  the  mean  of  each  data  set  had  to  be  accepteo  as  being 
equal  to  zero,  this  can  prove  to  be  misleading.  This 
result  should  be  easily  expected  considering  the  larte 
standard  deviations  in  the  data.  The  hypothesis  test  which 
contains  extremely  important  results,  from  the  tacticians 
viewpoint,  is  the  test  of  Section  lilt.  Remembering,  that 
standard  deviation  is  a  measure  of  the  precision  of  the 


prediction,  it  can  be  stated  quite  assuredly  that  the 
inversion  height  cannot  he  predicted  accurately,  either  by 
the  model  cr  persistence,  to  the  degree  required  for 
tactical  employment.  In  all  cases,  the  accuracy  of 
persistence  was  tetter  than  the  model,  but  still  r.ct  good 
enough  for  tactical  use. 

The  hypothesis  tests  concerning  the  equality  of  the 
means  revealed  that  there  is  no  significant  difference 
between  the  populations.  Once  again,  this  can  be 
misleading  due  to  the  large  standard  deviations 
enc  ountered . 

C.  PREDICTION  THROE  VS.  VERIFICATION  TIME 

Figure  6  clearly  indicates  the  model's  predictive 
ability  decreases  significantly  the  further  it  progresses 
into  the  24  hour  period.  Figure  7  shows  much  less  scatter 
associated  with,  the  persistence  predictions.  These  plots 
seem  to  indicate  that  the  model  was  expecting/predicting  a 
much  greater  rhange  in  Zi  than  actually  occurred. 

T.  PREDICTION  ERROR  VS.  MONTH  OF  OBSERVATION 

Figure  6  indicates  the  model's  predictive  ability  is 
significantly  reduced  during  the  summer  months.  There  is 
much  more  scatter  associated  with  the  summer  months  than 
others.  The  month  of  observation  had  no  discernible  affect 
on  persistence  as  seen  in  Figure  9. 


7.  SKY  CONDITIONS 

Initial  sky  conditions  had  no  apparent  affect  cn  the 
nod el  or  persistence. 

F .  GFNEPAL  CONCLUSIONS 

The  results  appear  to  point  toward  a  conclusion  that 
the  MA3L  model  and  persistence  are  not  useable,  from  a 
tactical  viewpoint.  More  testing  and  statistical  work 
should  be  done  before  a  judgement  can  be  made  about  the 
model  physics.  The  next  chapter  deals  with  approaches 
which  should  be  taken  in  an  effort  to  make  a  mere  valid 


evaluation  of  the  MABL  model 


VI.  FECO^ENEATICNS 


Proper  experiment al  planning  car.  give  a  reasonable 
assurance  that  the  results  of  an  experiment  will  provide 
clear-cut  answers  to  questions  under  investigation. 

The  examination  in  this  thesis  made  the  major 
assumption  that  the  radiosonde  data  was  "perfect"  and  that 
any  results  would  indicate  the  precision  of  the  models' 
predictive  ability.  However,  it  is  widely  known,  and 
accepted,  that  radiosonde  data  are  not  "perfect"  because 
the  radiosone  is  not  sensitive  enough  to  abrupt  changes  ir. 


atmospheric 

determine 


^^iitions  as  it  ascends, 
if  wie  error  contributed  by  the 


significant  and  if  the  radiosonde 


In  order  to 
radiosonde  is 
data  differs 


significantly  from  an  alternate  method  of  collecting  the 
same  type  of  data,  an  Analysis  of  Variance  ( ANCVA )  study 
should  be  done  involving  three  separate  comparisons.  A 
complete  discussion  of  this  method  can  be  found  in  Chapter 
12  of  Fef.  3.  The  mathematics  of  this  method  will  not  he 
discussed  here,  but  suggested  procedures  for  conducting  the 
necessary  data  collection  will  be  elaborated  upon. 

The  three  ANOVA  comparisons  recommended  are  delineated 
in  the  following  subsections. 


A.  ftA- IOSONTE  VS.  RADIOSONDE 


To  determine  if  the  data  received  by  ’’adioscndes  is 
consistent,  a  test  should  be  conducted  in  which  two 
radiosondes  are  launched  within  minutes  of  each  other  at 
each  observation  time .  Insuring  that  both  radiosondes  are 
identically  calibrated  and  launched  from  a  weather  deck 
where  there  will  be  little  possibility  of  collision  with 
other  structures  during  ascent,  the  first  radiosonde  should 
he  launched  and  then  the  second  launched  a  few  minutes 
later.  Dual  launches  should  be  made  at  each  measurement 
period.  The  number  of  dual  launches  made  should  he  limited 
to  between  5  and  12  because  of  cost  and  a  large  number  of 


’aunches 


will 


consi stency  in  da ta  . 


not  he  necessary  to  chpck  for 
It  is  imperative  that  the  A  NOV  A  tests 


he  carried  out  even  if  by  otserva  ion  the  data  appears 


consistent.  This  data  can  then  be  used  to  evaluate  the 


degree  o*  consistency  in  received  radiosonde  data. 
Regardless  whether  the  tests  determine  if  the  data  from  two 
different  radiosondes  is  consistent  or  net  consistent. 


further  tests  will  be  necessary. 


E.  RFFRACTOvETTR  VS.  RFFRACTOKETER 

The  AN/AMF-3  (XAN-5)  Airborne  Microwave  Refractometer 
(AMR)  is  currently  used  In  some  Navy  aircraft  to  obtain  the 
same  type  of  information  which  radiosondes  provide. 


3  i:iu  It  ateo’JS  d  at  a'  should  he  collected  by  two  aircraft  to 
determine  the  decree  of  consistency  in  ref  re 't ere  ter 
measurements.  This  test  will  require  the  same  aircraft  to 
carry  two  different  ref  Tactometers  or  two  aircraft  ea-.o 
carrying  an  Av?.  and  launched  within  minutes  of  each  'other 
and  flying  the  same  pattern.  As  a  result  of  the  A'-'fVA  test 
of  this  data,  if  the  consistency  is  goed  then  the 
F.ef ractometer  can  he  used  in  a  test  comparing  i*s'  data  to 
the  radiosonde. 

C.  FAEIOSONOE  VS.  REFRACTOKETE?. 

In  order  to  determine  if  there  is  any  significant 
difference  between  radiosonde  data  and  ref ract omet er  data, 
data  should  he  collected  within  70-6?  minutes  of  each 
other.  If  there  is  significant  difference  in  the  data  then 
careful  consideration  should  he  taken  to  determine  which, 
if  either,  method  provides  the  more  accurate  data. 

Once  these  tests  have  been  done  and  the  statistical 
results  determined,  then  much  clearer  statements  can  he 
made  about  the  model's  predictive  ability  because  more  will 
be  known  about  the  quality  of  the  input  data. 

T\  GENERAL  FECOvPEN DAT  IONS 

The  other  parameters  which  are  predicted  by  the  model 
Should  undergo  similar  statistical  analysis  as  that  which 
has  been  done  on  the  parameter  Zi.  This  analysis  is 


necessary  so  that  a  clearer  understanding  cf  the  model's 
predictive  abilities  is  in  own. 

Some  caution  and  consideration  should  he  given  to  the 
results  of  the  three  tests  suggested  above,  because  the 
results  for  any  or  all  of  the  tests  for  one  geographical 
location  and/or  tine  of  year  of  the  test  may  agree  or 
disagree  does  not  necessarily  mean  the  results  will  or  will 
not  agree  for  another  location  or  different  tire. 

*jc  attempt  should  be  made  to  alter  the  ^AEL  model 
until  the  A  NOV  A  tests  are  done.  A  qualitative  and 
quantitative  assessment  of  the  input  data  is  imperative 
before  a  fair  evaluation  of  the  MA?L  model  can  be  made. 


APPENDIX  A 


FP-S5  PROGRAMS 

The  programs  used  in  this  thesis  are  listed  on  the 
following  pages.  Each  program  can  he  used  on  an  FP-65  by 
first  starting  the  machine  and  placing  a  disk  with  the 
programs  on  it  into  the  disk  drive.  Next  type,  VASS 
STORAGE  IS  ":D?03",  then  hit  the  END  LINE  key,  then  type 
LOAD  ’’program  name"  and  then  the  ENT  LINE  key.  The  END  LINE 
key  must  he  used  after  an  entry.  Once  the  program  is 
loaded,  hit  the  RUN  key  and  the  program  will  prompt  for 
inputs . 

The  program  "EISTR"  will  first  prompt  for  a  data 
filename.  Upon  entering  the  filename,  the  next  prompt  will 
he  to  enter  the  number  of  data  points.  The  data  can  he 
entered  in  any  order.  Once  the  last  point  has  been  entered 
the  program  ranks  the  data  in  ascending  order,  computes  the 
mean,  standard  deviation,  maximum  and  minimum  values  of  the 
data,  and  the  range.  It  then  prints  the  data  plus  the 
results  of  the  computations. 

The  program  "PIN"  is  used  to  divide  the  data  into  tins 
of  equal  width  in  order  to  make  a  histogram  plot.  The 
program  prompts  for  the  data  filename,  then  the 
distribution  filename.  The  data  and  distribution  filenames 


must  diffe 


The  next  prompt  is  for  the  number  of 


This  value 
data  points 
there  are 


is  any  integer  value  less 
The  number  of  tins  chosen 
at  least  two  data  points 


h ins  . 

than  the  nurher  of 
should  he  such  that 
per  bin  whenever 


poss i ble . 

Upon  completion  of  the  "BIN'"  program,  "TPLCT"  can  he 
loaded  for  use  in  making  a  histogram  plot.  The  first 
prompt  is  for  the  data  filename,  however,  enter  the 
distribution  filename  for  which  a  plot  is  desired.  A 
prompt  for  the  title  and  subtitle  of  the  plot  are  displayed 
next.  If  either  or  both  of  these  are  not  required,  hit  the 
END  LINE  key  and  the  next  prompt  will  be  displayed.  The 
next  entries  are  X‘d  I N ,  XUAX  ,  YMI N ,  YMAX.  These  are  the 
maximum  and  minumum  values  for  the  horizontal  and  vertical 
axis,  respectively.  The  four  values  should  be  entered  or. 
the  same  line  separated  by  commas  before  keying  ENT  LINF 
The  next  prompt  is  X  AXIS  LOG  SCALE  LNTEP  Y  or  N.  Upon 
answering  yes  or  no,  the  next  prompt  is 
XINTV.YINTV, XINTC ,YINTC .  These  values  are  the  x  end  y 
intervals  for  labeling  of  the  x  and  y  axis  and  the  x  ar.d  y 
axis  intercepts.  The  x  and  y  axis  titles  are  requested 
next.  Upon  completing  these  entries,  a  plot  is  made  or.  the 
screen.  When  the  plot  is  done  hit  the  CONTINUE  key  and  the 
number  c  ffor  a  hard  copy)  and  it  will  plot  on  the  plotter. 
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eta  is  required  to  be  plotted.  The  first  prompt 
he  nurber  of  paired  date  sets.  The  next  query  is 
ate  fileoa~e.  The  next  pronpt  will  be  for  the  pei 
ets.  Enter  the  x  coordinate  first  followed  by- 
hen  the  v  coordinate  followed  by  the  EM  LINE  xey 
cmpletion  of  entering  the  data  sets,  "TPLOT"  can 


o  makp  a  scatter  plot. 


V  »T 


V  •  «- 


I 


» 


fa 


h 


♦ 


10  MASS  STORAGE  fS  *,0700* 

20  DISP  *  OAT  A  FILENAME* 

30  INPUT  D* 

40  0  I M  D (250  > 

50  DISP  *  NQ .  OF  DATA  POINTS* 
60  INPUT  0(0) 

70  S-(0(0)»1)  DIV  8 « 1 
80  B-8«8 


90  CREATE  Dt.16,0 

100  PRINT  -data  filename 

110  FOR  1 • l  TO  D ( 0 ) 

120  DISP  ’DATA  PT  * , I 
130  input  d(d 
140  NEXT  I 


150  FOR  J-l  TO  D ( 0) 
160  M-10000 


■>D* 


170  FOR  I  - J  TO  D ( 0 ) 

180  IF  OCIXM  THEN  M-D(I)  »  I- 
190  NEXT  I 

200  G-DCJ)  §  0 ( J) »M  9  D(L)-G 
210  NEXT  J 
220  S-0 
230  S2-0 


*  t u  r  UK 


»  u  u  v  u 


250  S-S‘0(I) 

260  PRINT  D ( I ) 
270  S2*S2»D< I ) *2 
280  NEXT  I 


I 


290  X-S/D(0) 

300  PRINT  ‘MEAN  VALUE 
310  Y-SQP(S2/0(0) -X“2) 
320  PRINT  *SOEV-  • ,y 
330  PRINT  * M I N  VALUE  • 
340  PRINT  *MAX  VALUE  • 


•»x 


i  D  ( 1 ) 
i  D  (D  (  0  )  ) 


350  PRINT 
) 


•RANGE  • 


v  j  |  vnr  m  ,  w 

370  FOR  1-0  TO  D( 
380  PRINT/  1  ,  D( 
390  NEXT  I 
400  ASSIGN/  1  TO 
410  END 


|D(D(0))-D(1 


J 

I 


1 


J 


.S 

J 


"j 


v 

I 
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T  P  L  D  T 


** 


10  REAL  Kl<8)  •  LOG  VALUES 
20  DIM  D ( 2500)  •  THE  DATA 
20  DIM  A 1 1 50  I , A 1  $ [ 5 0 ) ,8)1251  ,  Cl  I 
25  J 

40  GCLEAR 

50  DISP  'ANSWER  PROMPTS  WITH  'Y* 
OR  *  N  *  • 

60  P-0 

70  PLOTTER  IS  1 
80  GCLEAR 

90  IE  P-8  THEN  190  '  HARO  COPY 
IOC  IP  P-1  THEN  190  »  REDRAW 

110  IE  P-2  THEM  160  !  CHG  TITLES 

120  IF  P-3  THEN  180  1  CHG  X/Y  AX 

ES 

130  GOSUB  670  '  READOATA 

140  IF  P-4  THEN  200  •  MORE  DATA 

SAME  PLOT 

150  IF  P-5  THEM  200  '  MULTIPLOT 

160  GOSUB  380  •  TITLES 

170  IF  P-2  THEN  190  •  CHG  TITLES 

180  GOSUB  440  •  X/Y  AXES 

190  GOSUB  770  '  DRAW  AXES 

200  GOSUB  1430  •  PLOT  OATA 

210  PAUSE 

220  •  SELECT  PLOT 

230  OISP  'CHOOSE  PARAMETER* 

240  DISP  '  1-REDRAW* 

250  DISP  *  2-CHG  TITLES' 

260  DISP  '  3- CHG  X/Y  AXES' 

270  DISP  '  4-MORE  DATA  SAME 

PLOT* 

280  DISP  *  5-MULTIPLOT' 

290  DISP  '  6-NEW  PLOT* 

300  DISP  '  7-END  PLOTTING' 

310  DISP  '  8-HARD  COPY' 

320  INPUT  P 
330  IF  P-4  THEM  90 
340  IF  P-8  THEN  PLOTTER  IS  705  • 
GOTO  80 


350  IF  P-1  THEN  PLOTTER  IS  705  • 
GOTO  80 

360  IF  P-7  THEN  1680  !  END 
370  GOTO  70 
380  »  TITLES 
390  DISP  'TITLE' 

400  INPUT  Al 

410  DISP  'SUBTITLE' 

420  INPUT  All 
430  RETURN 
440  •  X/Y  AXES 

450  DISP  'INPUT  XMI N , XMAX , YMI N , Y 
MAX' 

460  INPUT  XI ,X2 ,Y1 ,Y2 

470  DISP  'X  AXIS  LOG  SCALE?'  .. 

480  DISP  'ENTER  Y  OR  N' 

490  INPUT  QO* 

500  IF  GOIO'Y'  THEN  560 
510  X3-0 
520  X4-X1 

530  DISP  'YINTV.YINTC' 

540  INPUT  Y3,Y4 
550  GOTO  660 

560  DISP  'INPUT  X  I NTV , Y I NTV , X I  NT 
C  VI HTC  • 

570  INPUT  X3.Y3.X4.Y4 

580  IF  P-0  THEN  620 

590  DISP  *CHG  AXIS  LABELS?  Y/N* 

600  INPUT  Oil 

610  IF  Qlt-'N'  THEN  660 

620  DISP  ' XAX I S  TITLE* 

630  INPUT  Bt 

640  DISP  ' YAX I S  TITLE' 

650  INPUT  C» 

660  RETURN 

670  !  READ  DATA 

680  DISP  'DATA  FILENAME' 

690  INPUT  FI 

700  ASSIGN#  1  TO  FI 

710  READ#  I  i  N9 

720  N9-2*N9 

730  FOR  1-1  TO  N9 


T  P  L  0  T  "  (cont.) 


740  READ#  1  i  DC!) 

750  NEXT  ] 

760  RETURN 

770  *  DRAW  AXES 

780  LOCATE  30,110,32,89 

790  CS12E  3 

800  •  DRAW  LINEAR** 

810  IF  QOtO  *Y*  THEH  GOSUB  1010 
820  •  OR  DRAW  SEN  I  LOG 
830  IF  Q0t-*Y*  THEN  GOSUB  1160 
840  >  PUT  LABELS 

850  PLOT  <X1*X2)/2,Y2*3«<Y2-Y1)/ 
26  -2 

860  LORG  6 
870  CSIZE  4 
880  LABEL  At 
890  LABEL  Alt 

900  PLOT  CX1«X23/2,Y1-(Y2-Y1)/B, 
-2 

910  LORG  4 
920  LABEL  Bt 
930  OEG 
940  LOIR  90 
950  LORG  6  ^ 

960  LORG  5  ■ 

970  PLOT  Xl-<X2*Xn/8,<Yl*Y2>/2, 
-2 

980  LABEL  Ct 
990  LOIR  0 
1000  RETURN 

1010  •  DRAW  LINEAR  AXIES 
1020  SCALE  X1,X2,Y1,Y2 
1030  AXES  X3.Y3.X4.Y4 
1040  LORG  5 

1050  FOR  X-Xl  TO  X2  STEP  X3 
1060  PLOT  X.Y1-CY2-Y13/55 
1070  LABEL  USING  *K*  ,  X 
1080  NEXT  X 
1090  LORG  8 

1100  FOR  Y-Yl  TO  Y2  STEP  Y3 
1110  PLOT  Xl-(X2-Xn/40,Y 
1120  LABEL  USING  *K*  |  Y 
1130  NEXT  Y 


1140  RETURN 

1150  •  ORAUI  LOGX-AXIS 

1160  LORG  6 

1170  IF  Lt»  "SET*  THEN  1210 
1180  Lt-*SET* 

1190  X1-LGTCX1) 

1200  X2-LGT(X2) 

1210  FOR  1-0  TO  8 
1220  X 1  C I > "LGT ( I  «2) 

1230  NEXT  I 

1240  SCALE  XI ,X2,Y1 ,Y2 

1250  MOVE  XI, Y1 

1260  DRAW  X2.Y1 

1270  FOR  I-Xl  TO  X2 

1280  MOVE  1 .Y1-CY2-Y13/55 

1290  LABEL  USING  *K*  ,  10*1 

1300  FOR  J-0  TO  8 

1310  K2-K1{J>»! 

1320  MOVE  X2.Y1 

1330  DRAW  K2,Y1*<Y2-Y1)/100 

1340  NEXT  J 

1350  NEXT  I 

1360  VAX  I S  XI, Y3 

1370  LORG  8 

1380  FOR  Y«Y1  TO  Y2  STEP  Y3 

1390  PLOT  X1-CX2-XD/40.Y 

1400  LABEL  USING  *K*  ,  Y 

1410  NEXT  V 

1420  RETURN 

1430  »  PLOT  DATA 

1440  DISP  'L1NETYPE  7* 

1450  DISP  *  1  SOLID* 

1460  DISP  *  2  END  PT  ONLY* 

1470  DISP  •  3  DOTS* 

1480  DISP  *  4  SDASH* 

1490  DISP  *  5  LDASH* 

1500  DISP  *  6  DASH  DOT* 

1510  DISP  *  7  LDASH  SDASH* 

1520  DISP  *  8  LDASH  SDASH  SDASH 

• 

1530  INPUT  L 
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T  ?  L  0  !  "  (cont.  ) 


1540  LINOTYPE  L 

1550  IE  QOtO'Y*  THEN  1500 

1560  MOVE  LGTCD(l)) ,0(2) 

1570  GOTO  1590 

1580  MOVE  DC1) ,D(2) 

1590  FOR  1-1  TO  M9  STEP  2 
1600  IF  OOIO'Y*  THEN  1630 
1610  PLOT  LGTCCC1)) ,DC1*1> 
1620  GOTO  1640 
1630  PLOT  D(!).CC!*1> 

1640  NEXT  1 
1650  PEMUP 
1660  PEMUP 
1670  PETURH 
1680  STOP 
1690  ENO 
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APPENDIX  2 

DELTA  M/DELTA  ?  VALUES 

The  Delta  v  and  Delta  ?  values  for  each  data  set  are 
listed  on  the  following  pages  in  the  output  format  of  the 
II P—35  "DISTF.”  program.  The  Delta  M  values  were  obtained  by 
subtracting  the  radiosonde  value  for  Zi  from  the  value 
predicted  by  the  model  at  that  time.  The  Delta  ?  values 
were  obtained  by  subtracting  the  current  radiosonde  value 
from  the  previous  value. 


DATA  FILENAME  MO  DELI 


-803 

-579 

-560 

-545 

-53S 

-422 

-421 

-323 

-272 

-252 

-234 

-234 

-214 

-211 

-211 

-1S4 

-179 

-170 

-124 

-112 

-97 

-59 

-84 

-83 

-70 

-50 

-28 

-15 

-7 

-5 

-2 

0 

12 

15 

16 

15 

35 

cc 

•w' 

73 

82 

81 

92 

111 

111 

132 

150 

159 

189 

221 

210 

225 

266 

285 

295 

321 

325 

351 

373 

394 

429 

446 

543 

549 

5*5 

552 

810 

£61 

1098 

1971 

MEAN  VAL'JE  =  80.ei2 

SDFV  =  410.026 

VIN  VALUE  =  -603 

MAI  VALUE  =  1971 

RANGE  =  2574 

DATA  FILENAME  MODELS 


R7Q 

53? 

583 

234 

211 

173 

112 

-84 

-53 

-7 

3 

16 

ss 

81 

111 

153 

212 

285 

325 

334 

543 

652 

256 


-552 

-422 

-272 

-224 

-211 

-170 

-9? 

-S3 

-26 

-5 

12 

15 

72 

92 

122 

189 

235 

295 

351 

429 

549 

810 


MEAN  VALUr  =  52.60 
SDEV  =  335.425 
VI N  VALUE  =  -579 
MAX  VALUE  =  1098 


DATA  FILENAME  M0DEL3 


-52= 

-401 

-170 

-124 

-64 

-63 

-50 

15 

73 

150 

MEAN  VALUE  =  -121.1 
S  DEV  =  196.091 
MIN  VALUE  =  -539 
.  VAX  VALUE  =  150 
HANOI  =663 


DATA  FILENAME  PERSISTS 

-464 

-214 

-155 

-75 

-30 

-12 

-4 

52 

159 

169 

MEAN  VALUE  =  -57.5 
SDEV  =  183.947 
MIN  VALUE  =  -464 
MAX  VALUE  =  159 


CATA  FILENAME  MCEEL5 


-234 

-214 

-211 

-70 

-25 

-2 

210 

255 

2S5 

44? 

543 

561 


MEAN  VALUE  =  155.31? 

TALUE"=  -234 
'1AX  VALUE  =  551 
RANGE  =1115 


rATA  FILENAME  PERSISTS 

-246 

-110 

-102 

-57 

-45 

93 

S3 

97 

115 

153 

190 

312 

MEAN  VALUE  =  40.917 
SDEV  =  146.575 
VI N  VALUE  =  -246 


DATA  FILENAME  MODELS 


-579 

-563 

-545 

-422 

-393 

-272 

-252 

-234 

-194 

-112 

-97 

-15 

-7 

16 

55 

111 

132 

159 

189 

235 

255 

351 

378 

354 

543 

555 

1098 

MEAN  VALUE  *  30.704 
SDEV  =  383.155 
MIN  VALUE  =  -579 
MAX  VALUE  =  1098 
RANGE  =  1577 


EAT A  7ILENANE  PERSISTS 


-762 
-577 
-418 
-355  ' 

-257 

-214 

-163 

-143 

-111 

-137 

-S3 

-42 

-25 

-17 

-a 

-3 

-2 

22 

32 

58 

59 
61 
92 

169 

177 

227 

ME/ N  VALUE  =  -88.519 
SEEV  =  223.998 
MIN  VALUE  =  -782 
MAX  VALUE  =  207 
RANGE  =  989 
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DATA  FILENAME  MODEL7 


-560 

-401 

-?S3 

-252 

-234 

-170 

-124 

15 

73 

80 

81 

92 

150 

159 

189 

210 

255 

325 

351 

394 

445 

543 

MEAN  VALUE  =  55.e*4 
SDEV  =  287.282 
MIN  VALUE  =  -560 
MAX  VALUE  =  543 


DATA  FILENAME  PERSIST? 


-7S2 
-214 
-214 
-153 
-1 55 
-114 
-111 
-110 
-75 
-57 
-30 
-23 
-17 
-3 
20 
54 
52 
15S 
159 
169 
177 
190 

MEAN  VALUE  =  -48.591 
SrCV  =  203.037 
MIN  VALUE  =  -782 
MAX  VALUE  *  190 
RANGE  =  972 


70 


DATA  FILENAME  PERSIST 


-57? 

-484 

-418 

-355 

-257 

-243 

-10? 

-102 

-89 

-45 

-42 

-29 

-25 

-12 

-9 

-5 

-4 

-2 

30 

32 

3? 

44 

59 

39 

61 

81 

93 

93 

9? 

115 

20? 

282 

213 

544 

565 

VEAN  VALUE 
SDEV  =  241 
MIN  VALUF  ’ 
MAX  VALUE  * 


=  0.0 
952 
:  -577 
555 


DATA  El LENAME  VCEEL9 


-194 

-112 

-83 

-70 

01 

xj 

111 

132 

285 

321 

429 

MEAN  VALUE  =  81.9 
S  DEV  =  193.228 
MIN  VALUE  =  -194. 
yAX  VALUE  =  429 
RANGF  =  823 


DATA.  EILENAME  PERSISTS 


-243 

-143 

-25 

40 

57 

71 

92 

99 

153 

158 


MEAN  VALUE  =  27.3 
SDEV  =  124.982 
VIN  VALUE  =  -246 
MAX  VALUF  =  156 
RANGE  =  416 


DATA  FILENAME  MODEL10 


-503 

-542 

-532 

-422 

-323 

-272 

-214 

-211 

-1S4 

-173 

-112 

-97 

-24 

-23 

-70 


0 

21 

159 

325 

579 

552 

610 

221 

1092 

MEAN  VALUE  =  22.92 
SDEV  =  450.575 
MIN  VALUE  =  -50 
MAX  VALUE  =  105 
F.ANG5  =  1701 


;ncj 


DATA  FILENAME  K0DEL11 


-560 

-401 

-252 

-234 

-234 

-211 

-170 

-124 

-CO 

-50 

-29 

-15 

-7 

-5 

12 

15 

15 

15 

35 

66 

73 

£0 

52 

111 

111 

132 

150 

189 

201 

210 

235 

265 

265 

295 

321 

351 

378 

394 

429 

446 

543 

549 

556 

1971 

MEAN 

VALUE  =  140.023 

STEV 

=  372.330 

MIN  VALUE  =  -550 

PAX  VALUE  *  1971 

FANG] 

:  =  2531 

EATA  FILENAME  PFRSIST11 


-7=2 

-356 

-249 

-214 

-214 

-153 

-155 

-114 

-lie 

-127 

-122 

-76 

-57 

-45 

-42 

-32 

-29 

-25 

-17 

-3 

_  R 

-4 

-3 

32 

32 

33 

4e 

44 

54 

5£ 

59 

61 

52 

71 

92 

93 

115 

169 

159 

177 

192 

227 

312 

1160 

NT  AN 

VALUE  =  7.356 

SEEV 

=  245.521 

y  I M 

VALUE  =  -752 

MAX 

VALUE  =  1152 

EANSE  =  1352 

APPENDIX  C 


PIcTOG?.A'iS  3E  PREDICTION  ERRORS 

The  following  pages  are  histograms  of  the  pr^d 
errors.  These  plots  shew  the  distribution  c-f  each 
set.  The  number  of  tins  were  chosen  such  that  there 
at  least  two  data  points  per  tin  whenever  possible. 


M0DEL2  PREDICTION  ERROR 
20  BINS 


PERSIST2  PREDICTION  ERROR 
20  BINS 


MODEM  PREDICTION  ERROR 
6  BINS 


PERSIST4  PREDICTION  ERROR 
6  BINS 


M0DEL5  PREDICTION  ERROR 
4  BINS 


PERSIST5  PREDICTION  ERROR 
4  BINS 


Noimndod  nib 


HD-A151  053 
UNCLASSIFIED 


STATISTICAL/TREND  ANALVSIS  OF  THE  NARINE  ATMOSPHERIC  2/2 

BOUNDARY  LAYER  MODEL (U>  NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY  CA  R  D  BISKING  SEP  84 

F/G  4/1  NL 


PERSIST11  PREDICTION  ERROR 
15  BINS 


REPRODUCED  AT  GOVERNMENT  EXPENSE 


a  s 


Nouvnndod  nig 
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appendix  r 


DELTA  M  -  DELTA  ?  VALUES 


The  following  page?  contain  the  data  sets  used  in  t 
paired  sarple  t  test.  These  values  were  attained 
subtracting  fror  Delta  v  the  corresponding  Delta  ?  value 
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r.aTA  FILENAME  PAIFFE1 


-423 

-2?? 

-225 

-130 

-84 

-54 

-22 

_e 

2? 

33 

40 

52 

88 

108 

155 

258 

322 

353 

402 

573 

1017 

MEAN  VALUE  =  106.594 
STEV  *  304.525 
MIN  VALUE  =  -493 
MAX  VALUE  =1156 
PA  MCE  =  1549 


1  0T3 


EAT  A  FILENAME  PA  IF. El 


-325 

-251 

-223 

-124 

-35 

-72 

-51 

-14 

1* 

25 

38 

44 

54 

37 

140 

150 

222 

2f7 

343 

356 

533 

784 


MEAN  VALUE  =  99.981 
SEEV  =  271. 39c 
MIN  VALUE  =  -4? 

MAX  VALUE  =  101 


DATA  FILENAME  PAIRZE4 

-403 

-272 

-235 

-225 

-165 

-132 

-85 

-51 

-14 

-5 

2^ 

2? 

33 

42 

52 

9? 

222 

272 

322 

349 

352 

365 

397 

565 

673 

791 

1017 

1155 

MEAN  VALUE  =  164.571 
STEV  =  377.412 
MIN  VALUE  =  -423 
MAX  VALUE  =  1156 


-AT A  FILENAME  PAIRED 


I 


-493 

-3?7 

-27? 

-225 

-155 

-133 

-139 

-95 

-?4 

-72 

-54 

-46 

-23 

-14 


17 

24 

27 


I 


28 

33 

43 


52 


64 

135 


,  143 

j  145 

•  153 

237 

323 

349 

673 

I  794 

1317 

1156 


MEAN  VALUE  =  94.171 

SE'EV  =  340.109 

I  MIN  VALUE  =  -493 

MAX  VALUE  =1156 
P.ANC-E  =  1549 


! 
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